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Abstract — In this paper, we propose a robust transceiver design 
for the A' -pair quasi-static MIMO interference channel. Each 
transmitter is equipped with M antennas, each receiver is 
equipped with N antennas, and the k tb transmitter sends Lk 
independent data streams to the desired receiver. In the literature, 
there exist a variety of theoretically promising transceiver designs 
for the interference channel such as interference alignment-based 
schemes, which have feasibility and practical limitations. In order 
to address practical system issues and requirements, we consider 
a transceiver design that enforces robustness against imperfect 
channel state information (CSI) as well as fair performance 
among the users in the interference channel. Specifically, we 
formulate the transceiver design as an optimization problem 
to maximize the worst-case signal-to-interference-plus-noise ratio 
among all users. We devise a low complexity iterative algorithm 
based on alternative optimization and semi-definite relaxation 
techniques. Numerical results verify the advantages of incor- 
porating into transceiver design for the interference channel 
important practical issues such as CSI uncertainty and fairness 
performance. 

Index Terms — Interference channel, robust transceiver, imper- 
fect CSI, preeoder design, decorrelator design, max-min fair, 
alternative optimization, semi-definite relaxation. 



I. Introduction 

In many wireless network scenarios, the channel is shared 
among multiple systems. The coexisting systems create mutual 
interference, which poses great challenges for communication 
systems design. Conventionally, interference is either treated 
as noise in the weak interference case U or canceled at the 
receiver in the strong interference case (2), (3|- In the past 
decade, various schemes are proposed to utilize multiple sig- 
naling dimensions for interference avoidance and mitigation. 
In particular, in the recent breakthrough work j4|, the authors 
show that the paradigm of interference alignment (IA) can 
be exploited to confine mutual interference to some lower 
dimensional subspace, so that desired signals can be trans- 
mitted on interference-free subspace. It is shown that this IA 
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scheme, if feasible, is optimal in the degree-of-freedom (DoF) 
sense. The results of |4| has triggered a number of extensions 
0, |6| and related works Q, (8). These IA-based schemes, 
albeit theoretically promising, have various limitations. First, 
IA-based schemes require ideal conditions to be feasible such 
as perfect channel state information (CSI) and very large 
dimensions on the signal space. For example, the conventional 
IA scheme |4| requires time or frequency extensions to have 
feasible solutions. For A' -pairs quasi-static MIMO interference 
channels where time / frequency extensions are not viable, the 
IA scheme [4] is only feasible for K < 3 (cf. J5)). Second, 
while IA-based schemes have promising DoF performance - 
which is an asymptotic performance measure for very high 
signal-to-noise ratio (SNR) - they are not optimal at medium 
SNR that correspond to practical applications. When designing 
practical communication systems for the interference channel, 
a number of technical issues shall be considered. Specifically, 
in practice only imperfect CSI is available and there are limited 
signaling dimensions. Moreover, it is important to ensure 
satisfactory performance among all the systems in the network. 

In this paper, we consider the problem of robust transceiver 
design for the AT-pair quasi-static MIMO interference chan- 
nel with fairness considerations. Specifically, 1) we apply 
robust design principles to provide resilience against CSI 
uncertainties; and 2) we formulate the transceiver design as 
a precoder-decorrelator optimization problem to maximize 
the worst-case signal-to-interference-plus-noise ratio (SINR) 
among all users in the interference channel. In the literature, 
precoder-decorrelator optimization for worst-case SINR are 
proposed for broadcast and point-to-point systems lfT0l - lfT3l . 
Specifically, in iflOl . Ifl3l the authors consider precoding 
design for the worst-case SINR in MISO broadcast channel, 
where it is shown that the preeoder optimization problem 
is always convex. In lfl2l the authors consider precoder- 
decorrelator design for the worst-case SINR MIMO broadcast 
channel using an iterative algorithm based on solving convex 
subproblems. On the other hand, in ifTTIl the authors consider 
a space-time coding scheme for the point-to-point channel 
with imperfect channel knowledge. However, these existing 
works cannot be extended to robust transceiver design for the 
MIMO interference channel, which presents the following key 
technical challenges. 

The Precoder-Decorrelator Optimization Problem is NP- 
Hard: The precoder-decorrelator optimization problem for the 
interference channel involves solving a separable homoge- 
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Fig. 1. System model. There are A" source-destination pairs where each 
source node is equipped with M antennas and each destination node is 
equipped with N antennas. The fc" 1 transmitter sends independent data 
streams to the desired receiver. 



neous quadratically constrained quadratic program (QCQP), 
which is NP-hard in general fl4l . |[T5l . One approach to 
facilitate solving this class of problems is to apply semidefinite 
relaxation (SDR) by relaxing rank constraints; this method 
was applied in precoding design for MISO broadcast channel 
OS), E) and for MISO multicast channel OH, ED- Although 
the resultant semidefinite program (SDP) may be solvable, the 
optimization in general does not always have the desired rank 
profile. 

Convergence of Alternative Optimization Algorithm: 

Our proposed solution is based on alternative optimization 
(AO). The method of AO was proposed in |20), ED for 
precoder and decorrelator optimization for multi-user MIMO 
broadcast channels. However, coupled with the rank con- 
strained SDP issues as well as the absence of uplink-downlink 
duality (as in the case of broadcast channels) El . l23l . 
establishing the convergence proof of the AO algorithm in the 
interference channel is non-trivial 1241 and traditional conver- 
gence proof l20l . EH cannot be applied to our situations. 

Notation: In the sequel, we adopt the following notations. 
R MxN , C MxN and Z MxN denote the set of real, complex 
and integer M x N matrices, respectively; R + denotes the 
set of positive real numbers; upper and lower case letters 
denote matrices and vectors, respectively; M N denotes the set 
of N x N Hermitian matrices; X >; denotes that X is a 
positive semi-definite matrix; (-) T and (•)' denote transpose 
and Hermitian transpose, respectively; rank(-) and Tr( ) denote 
matrix rank and trace, respectively; [X]( a ^ denotes the (a, b) th 
element of X; || • || denotes the Frobenius norm; !(■) denotes 
the indicate function; JC denotes the index set {1, . . . , K} and 
Ck denotes the index set {1, . . . ,£&}; 0^ denotes an N x 1 
vector of zeros and Ijv denotes an N x N identity matrix; 
E[-] denotes expectation; and CAT(/x, <&) denotes complex 
Gaussian distribution with mean fi and covariance matrix 3>. 

II. System Model and Review of Prior Works 

A. System Model 

We consider a MIMO interference channel consisting of K 
source-destination pairs where each source node is equipped 
with M antennas and each destination node is equipped with 



N antennas as shown in Fig. Q] For ease of exposition, 
we focus on the fc th user referring to source node Sk and 
destination node Dk', nevertheless, the same model applies to 
all other source-destination pairs. Specifically, Sk transmits 



to Dk, which performs 



Lk data streams s^ fe ' = [s^ 
linear detection. The received signal of Dk is interfered by the 
transmitted signals of all other users. To mitigate the impact of 
mutual interference, prior to transmission Sk precodes the data 
streams using the precoder matrix = [v^ . . ] <G 
C Mxifc and Dk decorrelates the received signal using the 
decorrelator matrix U (fc) = [u[ k) . . .u^] G C NxLk . It follows 
that the transmitted signal of Sk is given by 

x w = vW s w = Ef= fe iv i (fe) Si (fc) , 

the received signal of Dk is given by 

. n (k) 



(1) 



£f =1 H^')x« 
H (k,k) x (k) + ^k =i H (M x a) +n (fc) , (2 ) 



interference 

and the decorrelator output of is given by 

S«= (U<*>)V fc ) 



(3) 



desired signals 

■ Ef=i (uW)tH^) VW S « +(Tj (fe > ) V fc ) . 



where H (fcj) G C NxM is the fading channel from Sj to D k 
and n( fe ) - CAf(0 N , N Q I N ) is the AWGN. As per ©-©, the 

(k) 

estimate of data stream s ; is given by 



leakage interference 
-<7VxM 



( u W)t H ^^)vJ fe) S | fe) + E^ =1 K (fc) )tH 

v v ' m^l 



(k,k) (k) (fe) 
* in 



desired signal 



inter-stream interference 

Ef=i ^(uf^H^vH 5 ^ +(u| fc VnW, (4) 



leakage interference 

where the severity of the inter-stream and leakage interference 
terms depend on the transceiver processing and CSI assump- 
tion. Considering practical systems, we make the following 
assumptions towards designing effective precoders and decor- 
relators. 

Assumption 1 (Transmit power constraint): We assume the 
data streams are independent and have unit power, i.e. 
]E[( s ( fc ))t s O)] — Furthermore, we assume the maximum 
transmit power of the fc th source node is Pk so the pre- 
coders shall satisfy the power constraint E^x^^x^'] = 

Y.tM k) H k) < ■ 

Assumption 2 (Fading model): We assume quasi-static fad- 
ing so the fading channels H^' 3 -' remain unchanged dur- 
ing a fading block. In addition, we assume rank(H < -' c ' : '* 1 ) = 
min(M,A r ). ■ 

Assumption 3 (CSI model): We assume perfect CSI is 
available at the receivers (i.e. perfect CSIR), and only imper- 
fect CSI is available at the transmitters (i.e. imperfect CSIT) 
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for designing the precoders and decorrelators. Specifically, we 
model channel estimates at the transmitters as 



H 



CM 



H 



(M 



- A( fcj '\Vj,fce fC, 



(5) 



where A< fc ^ is the CSI error (TO), 03], ED. Specifically, 
we assume ||A^ fcj ^|| 2 < e, which implies that the actual 
channel belongs to a spherical uncertainty region cen- 



~(.k,j) 
tered at H 

we denote U = {H (fc ' j) }f fc=1 = {H 



{H (fcj) }f 



with radius e. For notational convenience, 
J, fAW)}f tl and 



j,k=V ~ 

Remark 1 (Interpretation of the CSI error model): The 

imperfect CSIT model (0 encapsulates the following 
scenarios. 

. Quantized CSI in FDD Systems fiM Section II-B]: For FDD 
systems, the transmitters are provided with quantized CSI via 
feedback. Using uniform quantizers, the quantization cells in 
the interior of the quantization region can be approximated 
by spherical regions of radius equal to the quantization step 
size. As a result, the imperfect CSIT model corresponds to 
quantized CSI obtained using a uniform vector quantizer with 
quantization step size %fe. 

• Estimated CSI in TDD Systems IHT] Section IV-A]: For 
TDD systems, the transmitters can estimate the channels 
from the sounding signals received in the reverse link. The 
imperfectness of the CSIT in this case comes from the 
estimation noise as well as delay. Using MMSE channel 

prediction, the CSI estimate H is unbiased, whereas the 
CSI error A^ k ' 3 ' is Gaussian distributed and independent 
from the CSI estimate H . As a result, A 1 * ' J> is a jointly 
Gaussian matrix and ||A'*' J ')|| 2 < e corresponds to "equal 
probability contour" on the probability space of &y*<i) . In 
other words, the probability of the event HAC^II 2 < e 
depends on £ only. Accordingly, we could find an e such 
that Pr[|| A( fc J) || 2 < e] = 0.99 (for example). ■ 

By Assumptions Q] to [3] the data stream estimate s?' in (0]i 
can be equivalently expressed as 



(fe,fe) 
(fch t /A(M0. 



sf>= (uf V(S +A (fc ' fc) )v, ( *M* ] 



(6) 



£m=l(Uj 

Ef=iE™'=i(u| fe) ) t (H (fcj) +A( fc ^)vH ) ^ ) + (u| fe) )tnW. 



The actual SINR of Sf ' at the /c* receiver is given by (O, 
whereby the instantaneous mutual information between data 
stream s[ and estimate can be expressed as 



Cf^{{vW}^jf =1 ,uW) 



(8) 



io g2 (i+ 7i (fc) (H,{{vwc =1 }f =1 ,ur ; ))- 



(.n-\ 



B. Review of Prominent Transceiver Designs for MIMO In- 
terference Channels 

In the following, we review the motivations and issues of 
prominent transceiver designs for MIMO interference channels 
in the literature. 



1) Interference Alignment in Quasi-Static MIMO Signal 
Space: In J4], f9) the authors exploited IA in quasi-static 
MIMO signal space for precoder-decorrelator design. Specif- 
ically, assuming perfect CSI, we could obtain precoders and 
decorrelators that confine the interference on each destination 
node to a lower dimension subspace, such that interference can 
be more effectively removed. Note that IA is only feasible with 
sufficiently large number of signaling dimensions. For the im- 
pair quasi-static N x M MIMO interference channel, IA could 
achieve a DoF of K """W^) f or JC < 3 but might not be 
feasible for K > 3. Moreover, IA is not optimal in general at 
medium SNR. For example, consider the data stream estimate 
sj fc ' ) in ©; suppose IA is feasible then 



(u| fc) )tfi 



,U) 



0, j k or I m, 



and the actual SINR of the I th data stream at fc th receiver is 
given by 



A" 



||( u «)t ( H<^) +A C*^)) v <*)||2 
7E1 



(9) 



+ Ej=iE^= 1 ll(»f') t A^v«)|| 2 +iV ||uW|p 

V &k ) 



As per (O, the presence of CSI error A^ fcj ' creates persistent 
residual interference. Even when the residual interference is 
negligible, i.e. 



K H 



|| (u (fc) ) t ( Q('--.fc) +A ( fc . fc ) )v (fc)||2 

JVol|u< fc) || 2 



the conventional IA scheme IH, (9| makes no attempt to 
optimize SINR performance. 

2) Interference Alignment in Real Fading Channels: In Q, 
151 the authors consider IA along the real line by creating 
fictitious signaling dimensions. Specifically, assuming perfect 
CSI, we could design the leakage interference terms at each 
destination node to have the same scaling factor (or pseudo 
direction), such that interference can be effectively removed. 
For example, consider the received signal in (O; for the 
purpose of illustration let M = N = 1 and H^-^ e K so 



H (kM) x (k) + j2f =1 H (k 



3)U3) 



- Ef-i(# (fc,fc) - 



ik)(k) 



3^k 



-&{kJ)) v 



leakage interference 

where (a) follows from (fl} and ©. To facilitate IA along the 

real line, the data streams shall belong to the set of integers 

(i.e. € Z) and we shall choose the precoders such that 
fj(k,j) v U) = ff(k, m ) v (m) for ■ ^ m ^ R is shown in gj^ 

that, if ideally CSI error is negligible (i.e. w 0), 

this scheme could theoretically achieve a DoF of K m+n ~ 
However, this scheme would require infinite SNR and cannot 
be implemented in practice. 
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\ K n {k) \ - 



l(u 



)yg> || 2 +Ef = i Em=i I (», (fc) ) f (H (fc ' J) + ACS )v$ 1 1 2+JV0 1 |u|' 



(7) 



3) Iterative Algorithms to Minimize Leakage Interfer- 
ence / Maximize SINR: In J6) the authors exploit 
uplink-downlink duality and propose iterative algorithms for 
precoder-decorrelator design. Specifically, the algorithms in 
J6] Algorithm 1], are established with the objective of 
sequentially minimizing the aggregate leakage interference 
induced by each data stream, whereas the algorithm in (6] 
Algorithm 2] is established with the objective to sequentially 
maximize the SINR of each data stream. Note that the afore- 
mentioned algorithms neglect the presence of CSI error, which 
could have significant performance impacts. Moreover, these 
algorithms neglect individual user performance and fairness. 
This is undesirable because for practical systems it is important 
to ensure all users have satisfactory performance. 

III. Problem Formulation: Robust Transceiver 
Design with Fairness Considerations 

In this section, we formulate a transceiver design for the 
A' -pair quasi-static MIMO interference channel that is robust 
against CSI uncertainties and with the objective of enforcing 
fairness among all users' data streams. Specifically, to provide 
the best resilience against CSI error, we adopt a worst-case 
design approach. On the other hand, the fairness aspect is 
motivated by the practical system consideration to ensure 
all users in the network can have satisfactory performance. 
As such, we formulate the precoder-decorrelator design with 
imperfect CSIT as an optimization problem to maximize the 
worst-case SINR among all users' data streams, subject to the 
maximum transmit power per source node. 

A. Optimization Problem 

The robust and fair transceiver optimization problem for 
the A'-pair N X M MIMO interference channel consists of 
the following components. 

• Optimization Variables: The optimization variables include 
the set of precoders {{Vm j-^LjJjLi an d tne set of decor- 
relators {{u^lmLij-jLi- These variables are adaptive with 

respect to imperfect CSIT ti = {H (fej) }^ fe=1 . 

• Optimization Objective: The optimization objective is to 
maximize, with imperfect CSIT, the minimum worst-case 
SINR amon£0 all users' data streams (perceived by the 
transmitter) given by (cf. © and Assumption [3} 

& ..a A. Tf'^a^jljf^.uf 1 ). do) 

viec k 

• Optimization Constraints: The optimization constraints 
are the maximum transmit power for each source node 

'Note that )10l is the worst-case SINR perceived by the transmitter based 

^ ^(k*j) j>- 

on imperfect CSIT H = {H ' }j fc=1 . We choose the worst-case SINR 
perceived by the transmitter in order to incorporate robustness against CSI 
error A^'-?). 



P±,...,Pk, which give the precoder power constraints 
Ef=\M fe TV; (fe) <P k ykeJC (cf. Assumption Q}. 
Accordingly, the optimization problem can be formally 
written as Problem 1. 

Problem 1: (Robust Max-Min Fair Precoder-Decorrelator 
Design): 

{{{(v^)*}^=i}f=i, {{(u«n™Li}f=i} = V(Pi, • . . , Pk) 

u u) eC «xiV*e£ fc 

s- 1. Et fc i(v| fc) )M fc) <n.,vfceA:.(iib) 

In (II lab , the worst-case SINR with imperfect CSIT is given 
by the following proposition. 

Proposition 1 (Worst-Case SINR with Imperfect CSIT): 
Given CSI estimates H = {H }jfc=i at tne transmitter 
with error ||A^ fcj ^|| 2 < e, the worst-case SINR of data stream 
estimate s| fc ' ) perceived by the transmitter can be expressed 
as (fL2l . 

Proof: Please refer to Appendix [A] ■ 
Using Proposition Q] and let P = min(Pi, . . . , Pk) and 
Pk = Pk/P, we can recast Problem V as 

{7^{{(v^)n^l}f=l,{{(uH ) r}^=l} J A :l} = 7'(^) 

min —7 (13a) 

v « gC Mxl 

□fee** 1 

7SR+ 

s. t. # ) (W,{{v^}^ 1 }jLi,uj fc) )>7,VlG-C fc ,VfcGX: ) (13b) 
£t\(v| fc) )M fe) < PkP,Vk e K. (13c) 

B. Properties of the Optimization Problem 

Note that it is not trivial to solve Problem V since it is 
non-convex and NP-hard in general as we elaborate below. 
In Section IIVI we shall propose a low complexity iterative 
algorithm for solving Problem V. 

1) Problem V is a non-convex problem: The minimum 
SINR constraints in d 1 3bb can be rearranged as 

(1 + -T)||( U p>)tlk (fc ' fc) V p^| 2 + l)e|| U p>|| 2 || V p^|| 2 

- 7 /Vol | U; (fc) || 2 -7 EjLl Em=l 1 1 (U^ )tfi (fc J) V« || 2 

-7^IK (fe) || 2 Ef=iE^=iHv^|| 2 >o, 

which are non-convex inequalities consisting of non-positive 
linear combinations of norms. Therefore, Problem V is a non- 
convex problem. 

2) Problem V is NP-hard in general: To illustrate that 
Problem V is NP-hard in general, we consider the inverse 
problem of jointly minimizing the transmit powers of all 
source nodes subject to a minimum SINR constraint for all 
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min_ 7 W(W,{{v«}^ =1 }f =1 ,uW) 



||A< fc -J) || 2 < 



(u W)ta"- fc >vW||'- e ||n|*>||'||v,W||' (12) 



z:5Lx i ic»J* J )*fi < *' : ° vS' i i a +=i i-J* 3 1 i a s:™*= a 1 1^ 1 1»- 1 1 cnf^j+fi 1 *** 5 ^* 1 i»J* 5 1 i a i k» w i i a +jv i i»f* 

ll \rLii\ y rn J- m=1 J-,- =1 ,Uj J. 



transmit 
power 



Problem V 



Target worst-case SINR 



precoders 
decorrelators 



Problem Q 



Optimize decorrelators 
with fixed precoders 
(Section IV-B) 

Problem Q u 



Optimize precoders 
with fixed decorrelators 
(Section IV-C) 

Problem Q v 



precoders 
decorrelators 



Fig. 2. Interrelationship among the optimization problems. 

users' data streams^. In Section IIVI we shall propose an 
algorithm for solving Problem V facilitated by solving the 
inverse problem^] that consists of the following components. 

• Optimization Variables: The optimization variables include 
the set of precoders {{Vm and the set of decorre- 
lators {{u«}^ =1 }f =1 . 

• Optimization Objective: The optimization objective is to 
minimize the required transmit power of all source nodes, by 
means of minimizing the precoder powers Ya^i^i ' 
Vfc G K. 

• Optimization Constraints: The optimization constraint is 
for all users' data streams to meet the prescribed minimum 
SINR 7 , i.e. 7l k \H,{{y$}%=i}f=iM k) ) > 7- 

Accordingly, the inverse problem can be formally written 
as Problem 2. 

Problem 2 (Power Minimization Precoder-Decorrelator De- 
sign): 

= 2(7) 

min j3 

v CO eC Sfxi 
/3GK+ 

s-t Ef=i(v ; (fc) )S (fc) < Pk p, Vfce/C, 



(14a) 



(14b) 



^(^{{vWl^xl^x^^T.V/er^Vfce/CXWc) 

Consider an instance of Problem Q with minimum SINR 
constraint 7, i.e. 



{PA{v\ih L rl=i}f=iA{^}2=i}f=i\ = Q(7), (is) 

and the required transmit power of the fc th source node is pk/3. 
It can be shown that 



{7 ! {{v^ ) }^ =1 }f =lJ {{S^C =1 }f =1 }=P( / 3) 



(16) 



2 Please refer to 1161 — 1181 , 1261 and references therein for discussions on 
the inverse relationship between max-min fair and minimum power precoder 
design problems for MISO broadcast and multicast channels. 

3 The inverse problem will be utilized in Section llV-CI 



so we can solve Problem Q to obtain a corresponding solution 
for Problem V, and vice-versa. Since Problem Q is NP-hard 
in general, Problem V is also NP-hard. Specifically, we define 
the special case of Problem Q with fixed decorrelators as 

Problem 3 (Power Minimization Precoder Design with Fixed 
Decorrelators): 

{e,{{<y$r}%=i}*=i} = 2v( 7 ,{{u«} I L3 =1 }f =1 ) 



min £ 

v « gC MXl 

<k 



(17a) 
(17b) 



7r^,{{^ ) }™ , = i}f=i,u| fc) )>7,V?e/: fc ,Vfce/C.(17c) 

Note that Problem Q v belongs to the class of separable 
homogenous QCQP, which is NP-hard in general lfT4l . lfT31 . 
This implies that Problem Q, which contains Problem Q v as 
special case, is also NP-hard in generafl 



IV. Low Complexity Iterative Solution 

In this section, we propose a low complexity iterative algo- 
rithm for solving the robust and fair transceiver optimization 
problem V . In particular, the proposed algorithm is facilitated 
by solving the inverse Problem Q, whereby we exploit the 
structure of Problem Q to apply effective optimization tech- 
niques. 

A. Overview of Algorithm 

The proposed algorithm for solving Problem V is facilitated 
by solving Problem Q as illustrated in Fig. |2j which is also 
detailed in Algorithm Q] Specifically, we iteratively refine 
the decorrelators and precoders to monotonically improve the 
minimum SINR. Each iteration consists of two stages: 
• (Steps 1-3 of Algorithm Q]) First, given the status quo mini- 
mum SINR 7 achieved with the fc th source node transmitting 
at power P^, we solve Problem Q to optimize the precoders 
and decorrelators for minimizing the transmit powers, i.e. 



{^.{{^^la^}^}^} = 6(7), (18) 

such that the minimum SINR 7 is achieved with the 
fc th source node transmitting at a reduced power of 



» (Steps 4-5 of Algorithm Q]) Second, we improve the min- 
imum SINR by up-scaling the transmit precoding powei0 
of the k lh user to the power constraint P/., i.e. v„ = 



4 Problem Q v will be utilized in Section IIV-CI 

5 We show in )3U that up-scaling the precoding powers improve the 
minimum SINR. 
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We repeat the iteration until the minimum SINR converges to 
a maximum. However, it is not trivial to solve the iteration 
step as per < fT8l > since Problem Q is NP-hard in general as 
shown in Section HII-B2I As such, we shall solve Problem Q 
based on alternative optimization between the decorrelators 
and the precoders, i.e. we present the algorithm for optimizing 
the decorrelators with fixed precoders in Section IIV-BI and 
introduce the algorithm for optimizing the precoders with fixed 
decorrelators in Section ITV-CI The top-level detail steps of the 
optimization algorithm is summarized below (Algorithm [TJ 
and illustrated in Fig [3] The convergence proof for AlgorithmQ] 
is provided in Appendix [D] 

Algorithm 1 (Top-Level Algorithm): 
Inputs: maximum transmit power for each source node 
P U ...,P K 

Outputs: precoders {{(Vm )*} m 3 = i}jLi ar >d decorrelators 

{{(u$)*}™U}f=i 

• Step 0: Initialize decorrelators {{S^}^L 1 }jL 1 and pre- 



coders {{v,„ } m J = i}j=ii where the transmit power for the 
j th source node is Em=x(^m = p r 

Repeat 

• Step 1: Optimize the decorrelators with fixed precoders (cf. 
Section HV^Bl i 

{{ u m'}m=l}jS=l = 2u({{ V m }m=l}j^l)- 

Update the candidate decorrelators (um)* = u^. 

• Step 2: Evaluate the minimum SINR 



mm7 i W (1W,{{vW}^ =1 }f =lJ Sr)=7. 

VfcG/C J 

Update the target SINR 7 = 7. 

> Step 3: Optimize the precoders with fixed decorrelators (cf. 
Section HV^Cl l 

{UivH^iljli} = Qv(7j{u^}^=i}f=i). 

> Step 4: Evaluate the required transmit power of each source 
node ft ^=E^=i(^ ) ) t vH ) - 

> Step 5: Evaluate the minimum SINR with up-scaled pre- 
coders 



vfce/c 
viec k 



7 and candidate precoders 



Update the target SIN R 7 

Until the minimum SINR 7 converges. 

Return precoders {{(\m )*}™=i}j=i an( ^ decorrelators 

{{(u«)*}^=Jf=i. 



B. Decorrelator Optimization with Fixed Precoders 

We define the decorrelator optimization problem with fixed 
precoders to maximize the minimum SINR among all users' 
data streams as 



^Initialize 1 



solve Problem Q„ 
(cf. Section IV B) 



Adjust target 
minimum SINR 



solve Problem Q, 
(cf. Section IV C) 



optimize decorrelators 
with fixed precoders 



optimize precoders 
with fixed decorrelators 



No 



'Minimum SINR' 
.converged? 



c 



Yes j 
Finish 



Fig. 3. Illustration of overall algorithm. 



Problem 4 (Maximum SINR Decorrelator Design with Fixed 
Precoders): 



{{(u^n^iijii = Q«({{v«}£ =1 }f =1 ) 

argmax min ^ k) (H, {{v«}^ =1 }f =1 , u ; (fc) ). (19) 



' ec Nxl 



viec k 



As per ( fl9l ), the worst-case SINR of data stream estimate 



s)^ only depends on decorrelator vS k . Therefore, we can 
independently optimize each decorrelator, i.e. 



(u w r = QW) ({{v a) } ^ i} ^ =l) 



argmax ^ k \u, {{vW}^ =1 }f =1 , u 



(20) 



and the optimal decorrelator is given by Theorem Q] 

Theorem 1 (Optimal Decorrelator with Fixed Precoders): 
Given precoders }m=i}$=i> tne optimal 



decorrelator for data stream estimate s\ is given by 

_ (F, w )-*(w} fc >" 



r(*0 



||(F< fc >)-2(wJ' 
A' tt^L 



where 



7=1 Z^m=l " 



£ Ej=l Em=l 



-e v 



hi 2 



|„0')||2t 
|Vm || In 



%'(v^) T (H 
H 



(fe '\«( v «)t(H^)t 



AT 



JVoI 



(Wj )* is the principle eigenvector of (Fj ) a E, (F, ) a , 

a n dEW=H (fc ' fe) v| fe )(v[ fe) )t(H (fe,fe V 
Proof: Please refer to Appendix [B] 



£ ||v, (fc) || 2 Uv. 



C. Precoder Optimization with Fixed Decorrelators 

In Section IIII-B2I we defined the precoder optimization 
problem with fixed decorrelators, Problem Q v (cf. ( 117a| i- 
d!7cl i). Since Problem Q v belongs to the class of separable 
homogenous QCQP, it is NP-hard in general. In the literature, 
some authors consider instances of this class of problems for 
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MISO broadcast channel that are always solvable (cf. IfTTl , 
||26l and references therein), whereas some authors consider 
problems for MISO multicast channel that are always NP- 
hard (cf. Ifl8l and references therein). For the interference 
channel model considered herein, we provide an algorithm 
for obtaining the optimal solution for Problem Q v . 

One effective approach for solving separable homogenous 
QCQP is to apply semidefinite relaxation (SDR) techniques. 
Let \{ k) = vpOpy. From O, the worst-case SINR of 
data stream estimate Sj can be expressed as (|2TI ). It follows 
that we can equivalently express the precoder optimization 
problem with fixed decorrelators as 

{SM{(V«)*}™U}f=i} = Qv(7,{{uH ) }™'=i}f = i) 

min S (22a) 

s-t Ef=iTr(v| fe) )<p fc H,VfceX;, (22b) 
r| fe) (H,{{ V« }% =1 }? =1 ,u| fc) ) > 7 ,VZ e £ fc ,Vfc e X(22c) 
V ; (fe) h 0,Vfc G £,VZ G £ fc , (22d) 
rank(V ; (fc) ) = l,Vfc G £,VZ G C k , (22e) 

where d22db and d22el i follow from the definition of Vj , 
( I22bl ) are power constraints, and ( I22cl ) are SINR constraints. 
Note that we could obtain the optimal precoder (Vm )* from 
the eigenvector of (V^)* corresponding to the only non-zero 
eigenvalue. 

Comparing between Problem Q v and Problem Q v , the SINR 
constraints of Problem Q v ( |22ct are convex inequalities, i.e. 

+( 7 - l)e| |u< fc) | | 2 Tr(V< fe } ) - 7 e| |u ; (fc) 1 1 2 £f =1 ^=1 Tr(vH } ) 

-7Ef= 1 E^=iTr((H (fcJ Vup ) (up ) )tH (M) V W) > 0, 

but Problem Q v is still a non-convex problem due to the 
rank constraints (I22eb . By means of SDR, we neglect the 
rank constraints and Problem Q v degenerates into an SDP 
that can be solved efficiently 11271 . In general, the resultant 
solution {S*, {{{Vm)*}m=i}jLi} could have arbitrary rank. 
If rank((V^ ) )*) = 1, Vm G Cj and Vj G /C, then constraints 
(122eb are intrinsically satisfied and {{(V^)*}^.^}!^ are 
optimal. The following theorem summarizes the optimality of 
the SDR solution in (l22at-d22eTi . 

Theorem 2 (Optimality of the SDR Solution): The SDR 
solution of Problem Q v will always give rank 1 solutions (i.e. 
rank((V£p)*) = 1) and hence, the SDR solution is optimal 
for Q v . 

Proof: Please refer to Appendix ICl ■ 

V. Simulation Results and Discussions 

In this section, we evaluate the proposed robust transceiver 
design via numerical simulations. In particular, we compare 
the performance of the proposed scheme against four baseline 
schemes: 

• Baseline 1 the conventional IA scheme J4); 

• Baseline 2 the SINR maximization scheme |6] Algorithm 2]; 



• Baseline 3 a naive max-min SINR scheme adopted from Ifl2l : 

• Baseline 4 a naive max-min SINR scheme adopted from (26). 
As discussed in Section IH-BI baselines 1 and 2 are theo- 
retically promising schemes for the interference channel but 
neglect important practical issues such as CSI uncertainty and 
fairness among users. On the other hand, baselines 3 and 4 
are adopted from existing max-min SINR schemes that are 
originally designed for the broadcast channel (i.e. there is only 
a single transmitter and multiple receivers). Without loss of 
generality, we assume independent and identically distributed 
(iid) Rayleigh fading channels, i.e. [H (fcj) ] (ajb) - CAf (0,1), 
Vj, k G JC, Va G [1,N], and Vfe G [1,M\. For the purpose 
of illustration, we consider the scenario where all users have 
the same power constraint P\ = . . . = Pk = P- In Fig. |4]to 
Fig. [6] we present simulation results for the average data rate^l 
versus SNB0 with different number of users and levels of CSI 
uncertainty. 

A. Fairness Performance 

In Fig.|4]and Fig. [5] we compare the average data rates of the 
proposed and baseline schemes. For the purpose of illustration, 
we consider the three-user 4x4 MIMO interference channel, 
where each user transmits L = 2 data streams and the pre- 
coders are designed with imperfect CSIT with e = {0.1, 0.15}, 
whereas the receivers have perfect CSIR. It can be observed 
that the proposed scheme achieves much higher average worst- 
case data rate per user than all the baseline schemes, and thus 
provides better minimum performance. For example, at CSI 
error e = 0.15, the proposed scheme has 5dB SNR gain over 
the SINR maximization algorithm (baseline 2) at providing 
a worst case data rate of 6 b/s/Hz and the conventional 
IA scheme (baseline 1) cannot provide worst-case data rate 
of 6 b/s/Hz. The superior performances of the proposed 
scheme is accountable to both the SDR approach as well 
as a suitably chosen utility function (optimizing the worst 
case performance). Specifically, the chosen utility function 
1) provide resilience against CSI uncertainties as well as 2) 
achieve fair performance among users. On the other hand, the 
SDR approach also contributes to obtaining a good solution 
for solving the optimization problem. 

B. Total Sum Data Rate Performance 

In Fig. [5] we compare the average total sum data rates of 
the proposed and baseline schemes for K = 3, N = M = 4, 
L = 2, and CSI error e = {0.1,0.15}. It can be observed 
that the proposed scheme not only achieves better worst-case 
data rate but also achieves higher total sum data rate than 
all the baseline schemes. In particular, due to the presence of 

6 The average data rate is defined as the average goodput (i.e. the bits/s/Hz 
successfully delivered to the receiver). Specifically, the goodput of data 
stream sf" is given by rp'l(r i '* : ' < c| fc '), where r\ ' = log 2 (l + 
7, (fe) ('H,{{v^ ) }^i}jii,u i (fc) )) is the scheduled data rate based on the 
SINR perceived with respect to imperfect CSIT H = {H ^'^ }j^ =1 , 

and Cf> = log 2 (l + 7i (fc) (W,{{vH , }^ =1 }ji 1 ,u i (fe) )) is the actual 
instantaneous mutual information. 

7 The SNR is defined as where No is the AWGN variance. 
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il k \n,{{v$}%=i}f=iA k) )= 



Tr((H ( ''- fc) )t„j fc '(„< fc ')tH< fc - fc )vj fc ))- £ || u ( fc )|PTr(vj fc) 



Ef =1 E^ 1 Tr((H <fc ' 3) )tu< fc )(u( fc »)tH< fc -' V U)) + e||u| fc »|| 2 i:f =1 E^ 1 TrK 



TrO^+JVolluHl* 



(21) 



r^£j{v«}^ = jf =1 ,u| fe >). 





6 8 10 12 14 16 18 20 22 24 26 28 

P/N (dB) 



Fig. 5. Average total sum data rate versus SNR. K = 3, N = M = 4, 
L = 2 and CSI error e = {0.1, 0.15}. 



algorithm (baseline 2) is designed assuming perfect CSI; its 
performance degrades rapidly for CSI error e > 0.02 and it 
can be observed that the achieved data rate could decrease 
with increasing SNR. On the other hand, the proposed scheme 
achieves a robust degradation with respect to CSI errors. 

VI. Conclusions 

In this paper, we proposed a robust transceiver design for 
the A' -pair quasi-static MIMO interference channel with fair- 
ness considerations. Specifically, we formulated the precoder- 
decorrelator design as an optimization problem to maximize 
the worst-case SINR among all users. We devised a low com- 
plexity iterative algorithm based on AO and SDR techniques. 
Numerical results verify the advantages of incorporating into 
transceiver design for the interference channel important prac- 
tical issues such as CSI uncertainty and fairness performance. 



CSI error, the total sum rate of the conventional IA scheme 
(baseline 1) does not scale linearly with the SNR anymore. 
Comparing Fig. [5] with Fig. |U it can be observed that the 
proposed scheme achieves the performance gain on fairness 
without sacrificing the total sum data rate. 

C. Robustness to CSI Errors 

In Fig. [6] we show the average worst-case data rates of 
the proposed and baseline schemes for different levels of 
CSI uncertainty. It can be observed that the proposed scheme 
always achieves higher average worst-case data rate than 
the baseline schemes. For example, the SINR maximization 



Appendix A 

Proof: Worst-Case SINR with Imperfect CSIT 

Given CSI estimates % = {H ' }j\. =1 at the trans- 
mitter, the worst-case SINR for each data stream estimate 
can be expressed as follows. Consider sj whose SINR 

l\ k \n, {{▼& ) }£ J = i}£ =1) u, (fc) ) is given by First, by the 

triangle inequality 

||( u p))t(H (fc ^ +A (M) v W||2 

> IKu^Vh^v^H 2 - ||(uj^A( fc >.>)v^|| 2 , 

||(u i (fc) )t(H (fcj) +A( fc J))vH ) || 2 

< IKu^Vh^v^H 2 + ||( u «)tA(M v «||2 ) 
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_ \\(u ( l k} )Hn ik ' k) +M k ^)y\ k) \\ 2 

~ E^Li ll(», (fc) ) t (H ( ' ! ' fc) +A(*.M)v^'||2+X;j c =1 2^ = i ||(u i (fc ')t(H (fc ' ;i) +A(fc.i))vH ) ll 2 +JVol|u 1 (fc) || 2 

|| (u W)tHt fc . fc ) v W|j2_|| (u (fc) ) t A ( fc . fc ) v W||2 

> || (u W ) tH"- fc )yW|P- e ||„W|P||vW|| 2 

^7« ( ' ,) («,{{^}S = i}^i.n l (,,) ). 



and so the SINR is lowered bounded as ( I23M . Second, with 
CSI error ||A^')|| 2 < e, 

||( u W )tA (M) v (j)||2 

= 1 ft((u 1 (fc) )tA( fc ^v^(v^)t(A^'))t Uj (fc) ) 



(a) 



< Tr(u i (fc) (uJ' c) )t)Tr(A( fc ^)v^ ) (vH ) )t(A( fcj ') t ) 



(6) 



< Tr(u i (fe) (u} fc) )t)Tr((A( fe ^)+A( fe ^)Tr(v^ ) (v^)t) 



= ||A( fc 'J')|| 2 



£ ||u, (fe) || 2 ||v«|| 2 , 



where (a) and (b) follow from the properties that Tr(AB) = 
Tr(BA) for A G C MxN and B G C NxM and Tr(CD) < 
Tr(C)Tr(D) for positive semi-definite C,D G C NxN . Thus, 
the worst-case SINR perceived by the transmitter can be 
expressed as (123 cl) . 

Appendix B 
Proof: Optimal Decorrelator with Fixed 
Precoders 

From dT2l . the worst-case SINR of data stream estimate sj fc ^ 
can be expressed as 

ll \rt, \\V m J- m=1 J- 3 - =1 ,Uj / — (u(' ! >)tF' fe) u' fs> ' 



where 



7=1 l^m=l ll V ' 



-£ V 



K 
3 

(k) 1 1 



0) 112 



Ijv-H 



(fe ' fe) v W(v«)t(H (fe ' fe) )t 



Ijv + A^oIjv, 



which is a Hermitian and positive definite matrix, and 



E 



(A') 



-H vj >(\\ >y(H y -e\\\\ >\\ Z I N , 



„Wl|2l 



which is a non-negative definite^ Hermitian matrix. Without 
loss of generality, let = c(Fj ) _ iw ; for arbitrary 

scaling factor c G C. We can equivalently expressed (l24l) as 



(»i t ')t(F' t ')4Ei t '(Ff')4,; t » 



(w}* J )tw}*' 



(25) 



(wj t) ) t QAQ t w | ' l ' ) 
(w< fc) )t«( fc > : 



8 If E (fe) 

is negative definite, then the CSI error e is too high. Without loss 
of generality, we assume e is sufficiently small. 



where QAQ^ denotes the eigen-decomposition of 

(Fp ) )-^ E / (fe) ( F | fc) )~^ 11 can be shown tha £l 42D is 
maximized with (Wj )* being the principal eigenvector 
of (Fj ) _ sEj (Fj fe ') — 3. In turn, the optimal unit norm 

decorrelator is given by (u| fc ^ ) 



||(FW)-i(wW).|| 

Appendix C 
Proof: Optimality of the SDR Solution for 
Problem Q v 

By using SDR, we solve the following SDP problem with 
complex-valued parameters: 



mm a 



(26a) 



s- t E^=iTr(V«)<^H,VjG/C 

Ef=iE^=iTr(A g$ V« ) > 6 J fc) , V/ G £ fc , Vfc G £ (26b) 
S > 0, (26c) 
V$ h 0, Vm G Cj , Vj G £, (26d) 



where A^'£%, G H M is given by 



L (Z,m)- 



(H ) T u ; '(a} ; ) T H — e||Uj || ] 
r (fc 'At..( fc V..( fc )Uw (fc J) i .1 i.,(*0 112 



if j = fc 
and m = l 



((H yup(uj J ) t H +e||u[ J || 2 I) otherwise 



and = 7A r o||u| fc ^|| 2 > 0. The corresponding dual problem 
is given by the following SDP: 



-axEf.iEfi^^ 

„(?) 



(27a) 



s. t. xWl-Ef =1 Efi A >: 0, VmG £ 3 - , VjG A(27b) 



i-Ef=i^ (i) >o, 

>0,Vfce/C,V/e£ fc , 
> 0,Vj G K, 

'Please refer to |28j Appendix E]. 



"As 



per 1291 Theorem 7.6.3] the principle eig' 



(27c) 
(27d) 
(27e) 

envalue of 



(F( fc ))-3Ej W (F, w ) _S is alwa Y s positive. 
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= x 



(k,j) 
(l,m) 



(29) 



+ Ef=i Ef=\ ((H (feJ Vu| fe) (u, (fc) )tH (fcj) + £ ||u| fc) || 2 l) I{fc + jM ^m} + y « ft£ ||u«|| 2 I 



rank M 



rank 1 



Note that (V$)* 7^ 0,Vj G A:,Vm G £,-, and from the 
complementary conditions for the primal and dual SDP: 



TrfZ^fV 



0, Vj G /C, Vm G 



(28) 



we can infer that ^ 0. Suppose that one of the optimal 



values {{(yP^j-^j-fcLi f° r me dual problem, say (y} l; ) 
0, then 



(IK* 



Z«=a; 



Wl+^yf'f-ftAfJlii ^ l&Z ^ 1} >- 0. 

fe=i 1=1 



It contradicts the fact ^ 0, and hence > 0, Vfc G 

K, V7 G £ fc . From d27bb and ([29), rank(Z^) > M — 1. On 
the other hand, from d28b , since Z$ ^Osorank(Z^) < M. 
It follows that rank(zH ) ) = M — 1. Moreover, due to ([28]) the 
optimal solution {{(V„ )*}^ l I =1 }^ 1 of primal problem d26l i 
must be of rank one. In other words, there will be zero duality 
gap between the primal non-convex problem Q v and the dual 
problem obtained by relaxing the rank constraint given by ( ITTI ). 

Appendix D 
Proof: Convergence of Algorithm Q] 

At the n th iteration of Algorithm Q] we denote 

Li 1 if 

3=1' 



the precoders as {{v^ [n]}^ J =1 }*L n , the decorrelators as 



{{u^' N}„ = i}jLi, the minimum SINR as j[n], and the 
transmit power scaling factor as j3[n). 

Upon initialization, we define the minimum SINR as 
7[0] = and start with arbitrary precoders {{vjjj [0]}m=i}j^=i> 
where the transmit power of the j th source node is 
Em=iK P])^m'M = Pj> anc ^ tne transmit power scaling 
factor is /3[0] = min (Pi, . . . , P K ). 

In the following, we show that each iteration of AlgorithmQ] 
increases the minimum SINR, i.e. 7[n] > j[n — 1], so 
Algorithm Q] must converge. 

In Step 1, given the precoders {{v,^[?i — l]}m=i}j^i' me 
decorrelators {{u,V? H}m = i}jLi are optimized to increase the 
minimum SINR, i.e. 



7= min (n, {{^ [n- 1] }^ =1 }f =1 , u\«> [n]) 
keic 



> min 7^ (H, [n- l]}^ =1 }f =1 , [n- 1]) 

k£lC 



(30) 



In Step 3, given the decorrelators {{u„ [n]}^ l J =1 }j^ 1 
and the minimum SINR constraint 7, the precoders 



{{Vm* [ n ]}m=i}j^i 316 optimized to jointly reduce the trans- 
mit powers of all nodes, i.e. the minimum SINR is unchanged 



7= min tT (U, [n]}^ =1 }f =1 , u\«> [n]) 

r r U) [ 



= min tT (H, {{v% [n-l]}^ =1 }f =1 , u\ K > [n]) 

kGlC 

whereas the transmit powers of all source nodes are reduced 



ft^N=E^=i(vH ) W) t v^[n] 



<E^Li(^[n-l])^[«-l] 
= Pj- 

In Step 5, the precod ers are up-sc aled to the power con- 
straint, i.e.Vm = \J Pjj (pj/3[n])\^m [n], where by defini- 
tion Pi/ pi = . . . = Pk/pk- As such, the minimum SINR is 
increased according to 

7 = min 7f > (H, {{v« [n]}£U}£i ^ [„]) 

k£K 



(J)r 



iec k ( 



|(ui fc 'w) t H (fc - fc) v," MII 2 - e ||u, (t) MII 2 ||v, ( '° ) MII 2 

Ef =1 E™U W)^* Nil 3 
fc) W II 2 Ef =1 Em=i I I'ff H ll 2 - £ l l u i <fc> M I l 2 l Iv| fc) M I 



V 

;( fc )/ 



- u n 



(*)[„])tH l ^* ) v< t > 



( fc '[n]|| 2 +JVo||u< fc »H|| 2 



/ 



<mm7r(^>{{V P ^/(^^N)^ ) N}^ =1 }f =1 ,ur ; W) 



t a(fe,fc> ()■) 



»]| 



|u, (fc) MII 2 ||v, (fc) NII 2 



kGlC 



7[nj 



E.^iE^IKurW)^^^!"]!! 2 \ 

el |u, (fc) [n] || 2 Ef =1 Em=i I 1^' N H 2 -e| |uf fc > [n] 1 1 2 | |v< fc) [n] 1 1 2 
. -||(u, <fc) [n])tH ( ' c - fc) v, (fc) [n]|| 2 + ^ |u| fc) [n]|| 2 , 

(31) 



It follows from ( f30b and ( |3l1 l that the minimum SINR increases 
with each iteration, i.e. j[n] > 7 > 7[n— 1], and Algorithm[T] 
must converge. 
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